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SCALAR CONSERVATION LAWS WITH MONOTONE 
PURE-JUMP MARKOV INITIAL CONDITIONS 


DAVID C. KASPAR AND FRAYDOUN REZAKHANLOU 


Abstract. In 2010 Menon and Srinivasan published a conjecture for the sta¬ 
tistical structure of solutions p to scalar conservation laws with certain Markov 
initial conditions, proposing a kinetic equation that should suffice to describe 
p(x , t) as a stochastic process in x with t fixed. In this article we verify an 
analogue of the conjecture for initial conditions which are bounded, monotone, 
and piecewise constant. Our argument uses a particle system representation of 
p(x , t) over 0 < x < L for L > 0, with a suitable random boundary condition 
at x = L. 


1. Introduction 


We are interested in the statistical properties of solutions to the Cauchy problem 
for the scalar conservation law 


p t =H(p) x in R x (0, oo) 
p = £ on K x {t = 0} 


( 1 . 1 ) 


where £ = £(x) is a stochastic process. In the special case 



( 1 . 2 ) 


( 0 ) is the well-known inviscid Burgers’ equation, which has often been considered 
with random initial data. Burgers himself, in his investigation of turbulence [8], 
was already moving in this direction, though it would be some time before problems 
in this area were addressed rigorously. 

Several papers, including have developed a theory showing a kind 

of integrability for the evolution of the law of p(-,t) for certain initial data. This 
article is concerned with pushing stochastic integrability results beyond the Burgers 
setting to general scalar conservation laws, an effort begun in earnest by Menon 
and Srinivasan [221125] . Before stating our main result, we first survey the major 
developments in this subject. 

This work is adapted from m- 
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1.1. Background. The Burgers case H(p) = —p 2 /2 has seen extensive interest. 
Recall that if 

u(x,t) = / p(x',t)dx' (1-3) 

•' — OC 

then u(x, t ) formally satisfies 

u t + = ® ( 1 - 4 ) 

and is determined by the Hopf-Lax formula m 

u{x, t) = inf |y £( 3 /)^'+ ^ |. (1.5) 


The backward Lagrangian y(x,t) is the rightmost minimizer in (11.51) . and it is 
standard [TBJ that this determines both the viscosity solution u(x,t) to (11.41) and 
the entropy solution to pt + pp x = 0, 

x-y(x,t ) 

P{x,t) = ---■ (1.6) 

Groeneboom, in a paper ms concerned not with PDE but concave majorant^ 
and isotonic estimators, determined the statistics of y(x,t) for the case where £(x) 
is white noise, which we understand to mean that the initial condition for u(x, t ) in 
(11.41) is a standard Brownian motion B(x). Among the results of [T5j is a density 
for 

sup{a: € R : B(x) — cx 2 is maximal} (1.7) 

expressed in terms of the Airy function. It has been observed (see e.g. [21]) that 
the Airy function seems to arise in a surprising number of seemingly unrelated 
stochastic problems, and it is notable that Burgers’ equation falls in this class. 

Several key papers appeared in the 1990s; we rely on the introduction of Cha- 
banol and Duchon’s 2004 paper m , which recounts some of this history, and discuss 
those portions most relevant for our present purposes. In 1992, Sinai m explicitly 
connected Burgers’ equation with white noise initial data to convex minorants of 
Brownian motion. Three years later Avellaneda and E [5] showed for the same 
initial data that the solution p(x, t ) is a Markov process in x for each fixed t > 0. 

Carraro and Duchon’s 1994 paper [9] defined a notion of statistical solution to 
Burgers’ equation. This is a time-indexed family of probability measures v(t, •) on a 
function space, and the definition of solution is stated in terms of the characteristic 
functional 

u(t,(j>) = J exp j p(x)(j)(x) dx^j i/(t,dp) (1.8) 

for test functions (j). Namely, u(t, <f>) must satisfy for each <f> a differential equa¬ 
tion obtained by formal differentiation under the assumption that p distributed 
according to v(t, dp) solves Burgers’ equation. This statistical solution approach 
was further developed in 1998 by the same authors m and by Chabanol and 
Duchon m It does have a drawback: given a (random) entropy solution p(x, t ) 
to the inviscid Burgers’ equation, the law of p(-,t) is a statistical solution, but it 
is not clear that a statistical solution yields a entropy solution, and at least one 
example is known [J] when these notions differ. Nonetheless, mm realized that it 


1 Straight forward manipulations of 11.51 ) relate this to a Legendre transform. 
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was natural to consider Levy process initial data, which set the stage for the next 
development. 

In 1998, Bertoin [3] proved a remarkable closure theorem for Levy initial data. 
We quote this, with adjustments to match our notation. 

Theorem 1.1 ( [4J Theorem 2]). Consider Burgers’ equation with initial data £(;r) 
which is a Levy process without positive jumps for x > 0, and £(x ) = 0 for x < 0. 
Assume that the expected value o/£(1) is positive, E£(l) > 0. Then, for each fixed 
t > 0, the backward Lagrangian y{x, t) has the property that 

y{x, t) — y(0, t) (1.9) 

is independent ofy(0,t ) and is in the parameter x a subordinator, i.e. a nondecreas¬ 
ing Levy process. Its distribution is the same as that of the first passage process 

x inf {2 > 0 : t£(z) + z > x}. (1.10) 

Further, denoting by ip(s) and Q(t,s) ( s > 0) the Laplace exponents of f(x) and 
y{x,t) - y(0, t), 

Eexp(s£(a;)) = exp(xip(s)) (1.11) 

Eexp [s(y(x,t) ~y(0,t))] = exp (x&(t,s)), (1.12) 

we have the functional identity 

s)) + Q(t, s) = s. (1.13) 


Remark. The requirement E£(l) > 0 can be relaxed, with minor modifications to 
the theorem, in light of the following elementary fact. Suppose that ^(ic) and 
f 2 (x) are two different initial conditions for Burgers’ equation, which are related by 
f 2 (x) = ^(a;) +cx. It is easily checked using on that the corresponding solutions 
p 1 (x,t) and p 2 (x,t) are related for t > 0 by 


p 2 {x,t) 


1 

1 + ct 


x t \ 
1 + ct’ 1 + ct) 


+ cx 


(1.14) 


This observation is found in a paper of Menon and Pego m, but (as this is elemen¬ 
tary) it may have been known previously. Using (11.141) we can adjust a statistical 
description for a case where E£(l) > 0 to cover the case of a Levy process with 
general mean drift. 


We find Theorem ll.ll remarkable for several reasons. First, in light of (11.61) , it fol¬ 
lows immediately that the solution p{x, t) — p{ 0, t) is for each fixed t a Levy process 
in the parameter x, and we have an example of an infinite-dimensional, nonlinear 
dynamical system (the PDE, Burgers’ equation) which preserves the independence 
and homogeneity properties of its random initial configuration. Second, the dis¬ 
tributional characterization of y[x, t ) is that of a first passage process, where the 
definition of y{x, t ) following (11.51) is that of a last passage process. Third, (11.131) 
can be used to show [23] that if if(t, q ) is the Laplace exponent of p(x , t ) — p(0, t), 
then 

+ V’V’g = 0 (1.15) 

for t > 0 and q € C with nonnegative real part. This shows for entropy solutions 
what had previously been observed by Carraro and Duchon for statistical solu¬ 
tions P33, namely that the Laplace exponent (11.121) evolves according to Burgers’ 
equation! 
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In 2007 Menon and Pego m used the Levy-Khintchine representation for the 
Laplace exponent (11.121) and observed that the evolution according to Burgers’ 
equation in (11.151) corresponds to a Smoluchowski coagulation equation ]EUE5] , with 
additive collision kernel, for the jump measure of the Levy process The 

jumps of y{-,t) correspond to shocks in the solution p(-,f). The relative velocity 
of successive shocks can be written as a sum of two functions, one depending on 
the positions of the shocks and the other proportional to the sum of the sizes of 
the jumps in ?/(•,£). Regarding the sizes of the jumps as the usual masses in the 
Smoluchowski equation, it is plausible that Smoluchowski equation with additive 
kernel should be relevant, and [El] provides the details that verify this. 

It is natural to wonder whether this evolution through Markov processes with 
simple statistical descriptions is a miracle m confined to the Burgers-Levy case, 
or an instance of a more general phenomenon. However, extending the results 
of Bertoin [1] beyond the Burgers case H(p) = — \p 2 remains a challenge. A 
different particular case, corresponding to L(q) = (— H)*(q ) = |g|, is a problem 
of determining Lipschitz minorants, and has been investigated by Abramson and 
Evans PQ. From the PDE perspective this is not as natural, since (—H)*(q) = \q\ 
corresponds to 

- H[p) = +ool(|p| > 1), (1.16) 

i.e. — H(p ) takes the value 0 on [—1, +1] and is equal to +oo elsewhere. So pQ, while 
very interesting from a stochastic processes perspective, has a specialized structure 
which is rather different from those cases we will consider. 

The biggest step toward understanding the problem for a wide class of H is 
found in a 2010 paper of Menon and Srinivasan Here it is shown that when 
the initial condition £ is a spectrally negativ ^ strong Markov process, the backward 
Lagrangian process ?/(•, i) and the solution p(-, t) remain Markov for fixed t > 0, the 
latter again being spectrally negative. The argument is adapted from that of H] 
and both U use the notion of splitting times (due to Getoor [Tl] ) to verify the 
Markov property according to its bare definition. In the Burgers-Levy case, the 
independence and homogeneity of the increments can be shown to survive, from 
which additional regularity is immediate using standard results about Levy pro¬ 
cesses [T5]. As [23] points out, without these properties it is not clear whether a 
Feller process initial condition leads to a Feller process in x at later times. Nonethe¬ 
less, [ 23 ] presents a very interesting conjecture for the evolution of the generator of 
which has a remarkably nice form and follows from multiple (nonrigorous, 
but persuasive) calculations. 

We now give a partial statement of this conjecture. The generator A of a sta¬ 
tionary, spectrally negative Feller process acts on test functions J £ C£°(M) by 

(AJ){y) = b(y)J'{y) + f ( J(z) - J(y)) f(y,dz ) (1.17) 

J — OO 

where b(y) characterizes the drift and f(y , ■) describes the law of the jumps. If we 
allow b and / to depend on f, we have a family of generators. The conjecture of |23j 
is that the evolution of the generator A for p(-,t) is given by the Lax equation 

A = [A, B\ = AB - BA (1.18) 


2 i.e. without positive jumps 
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for B which acts on test functions J by 

{BJ)(y) =-H'(y)b(t,y)J'(y)~ f - J(y)) f(t,y,dz). (1.19) 

J-oo V ~ z 

An equivalent form of the conjecture (11.181) involves a kinetic equation for /. The 
key result in the present article verifies that this kinetic equation holds in the special 
case we will consider. Before we state this, let us establish our working notation. 


1.2. Notation. Here we collect some of the various notation used later in the 
article. 

Write H[pi,... ,p n ] for the n th divided difference of H through pi,...,p n . For 
each n this function is symmetric in its arguments, and given by the following 
formulas in the cases n = 2 and n = 3 where Pi,P 2 ,P 3 are distinct: 


H\pi,P2] 


H(p2 ) - H(p 1 ) 
P2 - Pi 


and H[p 1: p 2 ,p 3 ) 


H[P2,P 3 \ -H\p 1 ,p 2 ] ^ 2Q) 

P 3 - Pi 


The definition for general n and standard properties can be found in several nu¬ 
merical analysis texts, including m- 

We write 5 X for the usual point mass assigning unit mass to x and zero elsewhere. 
For C > 0 and n a positive integer, write 

= {(ai,..., a n ) G R n : 0 < aq < • • • < a n < C) (1.21) 


and for the closure of this set in ffi”. We write dA% for the boundary of the 
simplex, and dtA^ for its various faces: 


d 0 A%-- 

= {(ai,.. 

. ,a„) S dA° : 

P 

II 




diA% = 

= {K-- 

., a„) S dA° 

■ Q>i — Q’i+l} 7 

* = !,-■ 

., n — 1 

(1.22) 

d n A c n - 

= {(«!,■• 

• ,a n ) S dA% : 

: a n = C}. 





Write M for the set of finite, regular (signed) measures on [0 , P], which is a Ba¬ 
nach space when equipped with the total variation norm || • || TV . Call its nonnegative 
subset M + . 

Let K, denote the set of bounded signed kernels from [0 ,P] to [0, P], i.e. the 
set of k : [0, P] i-A A4 which are measurable when [0, P] is endowed with its Borel 
cr-algebra and A4 is endowed with the er-algebra generated by evaluation on Borel 
subsets of [0, P], and which satisfy 

llfcll = sup{||fc(p_, -)I|tv : P- e [0, P]} < oo. (1.23) 

Observe that /C is a Banach space; completeness holds since a Cauchy sequence k n 
has k n (p _, •) Cauchy in total variation for each p_, and we obtain a pointwise limit 
A(p_, -)- Measurability of £:(•, A) for each Borel A C [0,P] then holds since this is 
a real-valued pointwise limit of k n (-,A). Write /C + for the subset of /C with range 
contained in Af+. 


1.3. Main result. In this section we provide a statistical description of solutions 
to the scalar conservation law when the initial condition is an increasing, pure-jump 
Markov process given by a rate kernel g. For this we require some assumptions on 
the rate kernel g and the Hamiltonian H. 
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Assumption 1.2. The initial condition £ = f(x) is a pure-jump Markov process 
starting at £(0) = 0 and evolving for x > 0 according to a rate kernel g(p-,dp+). 
We assume that for some constant P > 0 the kernel g is supported on 

{(p-,p + ):0<p.<p + <P} (1.24) 

and has total rate which is constant in p_: 

A = J g{p~,dp + ) (1.25) 

for all 0 < p- < P. In particular, g £ /C+. 

Assumption 1.3. The Hamiltonian function H : [0, P] — > R. is smooth, convex, 
has nonnegative right-derivative at p = 0 and noninfinite left-derivative at p = P. 


We provide for x > 0 a statistical description consisting of a one-dinrensional 
marginal at x = 0 and a rate kernel generating the rest of the path. The evolution 
of the rate kernel is given by the following kinetic equation, and the evolution of 
the marginal will be described in terms of the solution to the kinetic equation. 


Definition 1.4. We say that a continuous mapping / : [0, oo) —> /C+ is a solution 
of the kinetic equation 

{ ft = 

' (1.26) 


f(0 1 p-,dp+) = g(p-,dp+), 


where 

C*f(t,p-,dp+) 


= J( H iP*iP+\ - H[p-, p*])f(t, p-,d,p*)f(t, p*,d,p+) 

J H\p+,P*]f{t, P+, dp*) - J H[p-,p t ]f(t,p-,dp*) 

provided that 9 i-a C K f(9, •, •) £ K, is Bochner-integrable and 

I C K f(9,;-)d9 

J S 

for all 0 < s < t < oo. 


(1.27) 


f(t,p-,dp+), 


(1.28) 


Definition 1.5. We say that a continuous mapping £ : [0, oo) —>• M+ is a solution 
of the marginal equation 

j it = C°£ 

\^(0,dp o ) = So(dpo), ^ ^ 

where 


C°£(t,dpo) = / H[p*,p 0 ]£(t,dp*)f(t,p*,dpo ) 


H[p 0 ,P*]f(t,Po,dp *) 
provided that 9 i-a C°£{9, •) £ A4 is Bochner-integrable and 

£(t, ■) - £(s, ■) = [ C°£(6, •) d9 

J S 


£(t,d Po ), (1.30) 


(1.31) 


for all 0 < s < t < oo. 
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Theorem 1.6. Suppose Assumptionsand [~Q1 hold. Then the kinetic and mar¬ 
ginal equations have unique solutions, in the sense of Definitions [72)1 andlTThl Fur¬ 
thermore, the total integrals are conserved: 

A = J f(t,p-,dp+) 

(1.32) 

O 

II 

i — i 

(1.33) 

for all t > 0 and all 0 < p- < P. Finally, if 

g{p~i [0, P-] u {R}) = 0 

(1.34) 


for all 0 < p- < P, then f(t, ■, •) has the same property for t > 0. 

The kernels described by the Theorem [TTG] are precisely what we need to describe 
the statistics of the solution p 1 which brings us to our main result: 

Theorem 1.7. When Assumptions \1.2\ and 1 1. 3\ hold, the solution p to 

(p t =H(p) x (x,t) e R x (0,oo) 

\ P = f (i,t)eix{t = 0} (1 ' 35) 

for each fixed t > 0 has x = 0 marginal given by (.{t, dpo) and for 0 < x < oo evolves 
according to rate kernel f(t,p-,dp+). 


Remark. Theorems [Oil and H771 establish rigorously the conjectured 125. evolution 
according to the Lax pair (11.181) , within the present hypotheses. See [251 Section 2.7] 
for a calculation showing the equivalence of the kinetic and Lax pair formulations, 
which simplifies considerably in the present case due to the absence of drift terms. 
The Lax pair and integrable systems approach (in the case of finitely many states, 
where the generator is a triangular matrix) have been further explored by Menon 
m and in a forthcoming work by Li EQ3- 


1.4. Organization. The remainder of this article is organized as follows. In Sec¬ 
tion [2] we show that Theorem o will follow from a similar statistical characteri¬ 
zation for the solution to the PDE over x £ [0, L\ with a time-dependent random 
boundary condition at x = L. The latter we can study using a sticky particle system 
whose dimension is random and unbounded, but almost surely finite. Elementary 
arguments are used to check that our candidate for the law matches the evolution 
of the random initial condition according to the dynamics. Next in Section [8] we 
show existence and uniqueness of the solutions to marginal and kinetic equations 
of Theorem II.(il which are needed to construct the candidate law. The concluding 
Section 0] indicates some desired extensions and similar questions for future work. 
To keep the main development concise, proofs of lemmas have been deferred to 
Appendix [A] 


2. A RANDOM PARTICLE SYSTEM 

As functions of x, the solutions p(x , t) we consider all have the form depicted in 
Figure [T] From a PDE perspective this situation is standard—a concatenation of 
Riemann problems for the scalar conservation law—and we can describe the solution 
completely in terms of a particle system. Each shock consists of some number of 
particles stuck together, and the particles move at constant velocities according to 
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the Rankine-Hugoniot condition except when they collide. The collisions are totally 
inelastic. 


P 

A 


P4 


P 3 


P2 

Pi 


PQ 


Xl 


X 2 X 3 


X4 


Figure 1. For each t > 0, the solution p(x,t) is a nondecreasing, 
pure-jump process in x. We will see that for any fixed L > 0, we 
have a.s. finitely many jumps for x £ [0, L] and that p(-,t) on this 
interval can be described by two (finite) nondecreasing sequences 
(xi, ■ ■ ■ ,Xn', po, ■ ■ ■ ,Pn )■ 


The utility of this particle description is lessened, however, by the fact that 
the dynamics are quite infinite dimensional for our pure-jump initial condition 
£(x), x £ [0, + 00 ). To have a simple description as motion at constant velocities, 
punctuated by occasional collisions, we might argue that on each fixed bounded 
interval of space we have finitely many collisions in a bounded interval of time, and 
piece together whatever statistical descriptions we might obtain for the solution on 
these various intervals. Inspired by those situations in statistical mechanics where 
boundary conditions become irrelevant in an infinite-volume limit, we pursue a 
different approach. We construct solutions to a problem on a bounded space interval 
[0, L], and choose a random boundary condition at x = L to obtain an exact match 
with the kinetic equations. The involved analysis will all pertain to the following 
result. 

Theorem 2.1. Suppose Assumptions HOI a,nd M.i-ft hold. For any fixed L > 0, 
consider the scalar conservation law 

fp t =H(p) x (x, t) £ (0, L) x (0, 00 ) 

J p = £ x £ [0, L\ x {t = 0} 

[ p = ( (x, t) £ {x = L} x (0, 00 ) 


( 2 . 1 ) 
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with initial condition f (restricted to [0 ,L\), open boundar ^ at x = 0, and random 
boundary ( at x = L. Suppose the process ( has C(0) = £(L) and evolves according 
to the time-dependent rate kernel H[p t p + \f(t, pdp + ) independently of f given £(L). 
Then for all t > 0 the law of p(-, t) is as follows: 

(i) the x = 0 marginal is £(t, dpo), and 

(ii) the rest of the path is a pure-jump process with rate kernel f(t 1 p-,dp+). 

To prove our main result we can send L —> oo, applying Theorem 12.11 on each 
[0, L], and use bounded speed of propagation to limit the respective influences of far 
away particles (unbounded system) or truncation with random boundary (bounded 
system). The argument is quite short. 

Proof of Theorem M.A Fix any t > 0. We will write p(x, t ) for the solution to (11.351) 
over the semi-infinite a;-interval [0,+oo) with initial data £(#)• We take the right- 
continuous version of the solution. For L to be specified shortly, write p(x, t) for 
the solution and f,(x) for the initial data corresponding to (12.11) . 

Fix any x\,...,x k £ [0,+oo), and let A' = maxjxi,...,^}. Choose L > 
X + tH'(P). We couple the bounded and unbounded systems by requiring 

£I[o,l]=£, (2.2) 

allowing the random boundary ( to evolve independently of f given f (L) 

Recall m that the scalar conservation law has finite speed of propagation. Our 
solutions are bounded in [0, P], so the speed is bounded by H'(P). Since p{x,0) 
and p(x, 0) are a.s. equal on [0, L], we have also p(x,t) = p{x,t) a.s. for x £ 
[0, L — tH'(P)] D [0, X ]. From this we deduce the distributional equality 

{p(x 1 ,t),...,p(x k ,t)) = (p(xi,t),...,p(x k ,t)). (2.3) 

By Theorem 12.11 the latter distribution is exactly that of a process started ac¬ 
cording to £(t,dpo), evolving according to rate kernel f(t,p-,dp+). This process 
is terminated deterministically at x = L, which does not alter finite-dimensional 
distributions prior to x = L. Since p(x , t) is right-continuous and has the correct 
finite-dimensional distributions, the result follows. □ 

2.1. The dynamics. Our work in the remainder of this section is to prove Theo¬ 
rem IQ We begin by describing precisely those particle dynamics which determine 
the solution to the PDE. 

Figure [T] illustrates a parametrization of a nondecreasing pure-jump process on 
[0, Too) with heights po. Pi, P2, ■ ■ ■ and jump locations x\, X2, • ■ ■■ Our sign restric¬ 
tion on the jumps excludes rarefaction waves, and we have only constant values 
separated by shocks. Going forward in time, the shocks move according to the 
Rankine-Hugoniot condition 

±i =-H[pi-i,pi\ (2.4) 

until they collide. We say that each shock consists initially of one particle moving 
at the velocity indicated above. Then result of a collision can be characterized in 
two equivalent ways: 

^Using Assumption 11.31 and >0, the shocks and characteristics only flow outward across 
x = 0. Any boundary condition we would assign, unless it involved negative values, would thus 
be irrelevant. 
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(i) At the first instant when Xi = Xj+i, the i th particle is annihilated, and the 
velocity of the (i + l) th particle changes from —H[pi, pi + 1 ] to 

Xi+i =-H[pi-i,p i+ i]. (2.5) 

In the cas^l where several consecutive particles collide with each other at 
the same instant, all but the rightmost particle is annihilated. Since we 
only seek a statistical description for x > 0, we annihilate the first particle 
when X\ = 0, replace po with pi, and relabel the other particles accordingly. 

(ii) When Xi-i < Xi = x^+i = •■• = Xj -1 = Xj < Xj+i, the particles i 
through j all move with common velocity —H[pi—±,pj]. Following Brenier 
et al 0E] we call these a sticky particle dynamics. We can additionally 
take the position x = 0 to be absorbing: X, = 0 whenever x* = 0. 

We adopt the first viewpoint for convenience, as this is compatible with [12] (see 
the text following Definition 12.31 below), but a suitable argument could be given for 
the second alternative as well. 

Definition 2.2. For L as in Theorem 12.11 the configuration space Q for the sticky 
particle dynamics is 

OO 

Q=[]Qn, Qn=A^X A^. (2.6) 

n=0 

A typical configuration is q = (x±,.... x n ; p o,..., p n ) £ Q n when n > 0, or q = 
(po) G Qo = {Po ■ 0<Po<P} when n = 0. 

Definition 2.3. Our notation for the particle dynamics is as follows: 

(i) For 0 < s < t and q £ Q, write 

4>U = 4>Vq (2.7) 

for the deterministic evolution from time s to t of the configuration q ac¬ 
cording to the annihilating particle dynamics for the PDE, without random 
entry dynamics at x = L. 

(ii) Given a configuration q = (xi,..., x„, po,..., p n ) and p+ > p n , write e p+ q 
for the configuration (xi,..., x n , L, po ,..., p n , p+). 

(iii) Write for the random evolution of the configuration according to de¬ 

terministic particle dynamics interrupted with random entries at x = L ac¬ 
cording to the boundary process ( of (12.11) , where the latter has been started 
at time s with value p n . In particular, if the jumps of £ between times s 
and t occur at times s < T\ < ■ ■ ■ < < t with values p n +i < ■ ■ ■ < p n +k , 

then 

= 4>^Pn +k 4> T rt_^P n+k -r ■ ■ ■ <t>ll e Pn + 10? 9- (2-8) 

Proposition 2.4. For any s, q, the process <I>‘g is strong Markov. 

This assertion follows after recognizing as a piecewise-deterministic Markov 
process described in some generality by Davis [12] , Namely, we augment the config¬ 
uration space Q to add the time parameter to each component Q ni and then have 
a deterministic flow according to the vector field 

(1; -H[p 0 , pi],..., -H[p n - 1, Pn]; 0,..., 0). 

^which will almost surely not occur for the randomness we consider 


(2.9) 
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With rate f H[p n , p+\f(t, p n , dp+) we jump to the indicated point in Q n + 1 , and 
upon hitting a boundary we transition to Qk for suitable k < n, annihilating 
particles in the manner described above. 


2.2. Checking the candidate measure. Our goal is to take q distributed ac¬ 
cording to the initial condition and exactly describe the law of q for each T > 0. 
Using the kinetic (11.271) and marginal equations (11.301) . we construct for each time 
t > 0 a candidate law p(t, dq) on Q as follows. Take N to be Poisson with rate XL, 
Xi, ...,Xpf uniform on A^, and p 0 ,... , p^ distributed on Aj^ +1 according to the 
marginal t and transitions / independently of the xp. 

OO 

e~ XL J2 s n(d N )p n (t,dq), (2.10) 

n—0 

where po(t,dq) = l(t,dpo) and 

n 

H n (t,dq) = 1 a l(xi, .. . ,x n )dxi ■ ■ ■ dx n e(t,dp 0 ) Y[f(t,Pj-i,d Pj ), (2.11) 

3=1 

for n > 0. When we have verified (11.321) . it will be immediate that the total mass 
of this is one. 

We decompose the mapping from configurations q to solutions (as functions of 
x over [0, L\) into the map from q to the measure 

n 

n(q,dx) = poSo + ^ipi-Pi-i)S Xi (q € Q n ) (2-12) 

i-1 


and integration over [0,cc]. We claim that when q is distributed according to 
p(0,dq), the law of the measure 7r(<hg<7, •) is identical to that of ir(q',-) where q' 
is distributed according to p,(t, dq'). 

We now describe the structure of the proof of Theorem l2.ll Fix some time T > 0 
and consider F(t, q) = E q) where G takes the form of a Laplace functional: 


G(q) = exp j J{x)Tr(q,dx)^ = exp ^-p o J(0) - 'Y'XPi ~ Pi-i)^(a 


for J > 0 a continuous function on [0, L\. We aim to show that 

j t J ®G($Tq)n{t,dq) = 0 

for 0 < t < T, from which it will follow that 


(2.13) 


(2.14) 


J EG($jfr)/i(0,dg)= J G{q)p(T,dq) 


(2.15) 


for all G of the form in (12.131) . Using the standard fact that Laplace functionals 
completely determine the law of any random measure m, this will suffice to show 
that law of q,dx) for q distributed as p(0,dq) is precisely the pushforward 
through 7r of p{T, dq), and obtain the result. 

We might verify (|2.14p by establishing regularity for F(t,q) = EG(4>^ q) in t 
and the x-components of q. We speculate that it should be possible to do so if 
J is smooth with J'(0) = 0 and we content ourselves to divide the configuration 
space into finitely many regions, each of which corresponds to a definite order 
of deterministic collisions in </>. On the other hand, the measures p(t, dq) enjoy 
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considerable regularity (uniformity, in fact) in x, and we pursue instead an argument 
along these lines. Here continuity of F(t, q) in t will suffice. 

Lemma 2.5. Let G take the form of (12.1311 . Then F(t , q) = E q) is a bounded 
function which is uniformly continuous in t uniformly in q: 

w{5) = sup{|F(f, q) - F(s,q)\ : s,t G [0,T], \t - s| < 6, q G Q} ->■ 0 (2.16) 

as 6 —)■ O-l - . 

To differentiate in (12.141) we will need to compare EG(d>;fg) and E G($fq) for 
0 < s < t < T. Our next observation is that when the t — s is small, we can 
separate the deterministic and stochastic portions of the dynamics over the time 
interval [s, t]. The idea is that in an interval of time [s, t], the probability of multiple 
particle entries is 0{{t — s) 2 ), and a single particle entry has probability 0(t — s ) 
which permits additional o(l) errors. 

Lemma 2.6. Let 0 < s < t < T and q G Q n . There exist a random variable 
t G (s,t) a.s., with law depending on s,t, and q only through p n , and a constant 
C\,h'(p) independent ofq,s,t so that 


F{s,q) - F(t, 4>\q) - ( t-s ) J[EF(t,e p+ (f> T s q) - F(t, q)\H[p n , p+]f(t, p n , dp+) 

< - s ) 2 + (t- s)w(t - s)]. (2.17) 

We proceed to an analysis of the deterministic portion of the flow, tflq, where 
q is distributed according to p,(t,dq). For this we introduce some notation. Given 
q G Q n , we write 


a(q) = min-jr > s : (f> r s q G (dA%) x A^ +1 }. (2.18) 

For each fixed n, we consider several subsets of the set of n-particle configurations 
Q n . In particular we need to separate those configurations which experience an 
exit at x = 0 or a collision in a time interval shorter than t — s. For i = 0,. .., n — 1, 
write: 

A = {q G Q n ■ cr{q) < t and (f a s {q) q G (<%A£) X A^ +1 } 

n—1 

U = Qn \ 1J - 4 ' (2.19) 

B = {q G Q n : x n > L - (t - s)H[p n -i, p n }} 

V = Q n \B. 

Figure [2] illustrates these sets in the case n = 2. We observe in particular that (jfJJ 
and V differ by a p, n {t, -)-null set, and that the terms associated with the boundary 
faces are related to configurations with one fewer particle. 

Lemma 2.7. For each positive integer n and with errors bounded by 


CL,n,\,H'(P)[(t ~ s ) + w {t ~ s )] 


(2.20) 
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we have the following approximations: 




F ^ t ^ H[p*, p 0 \£(t, dp*)f(t, p*,dp 0 ) 


£{t, dp 0 ) 


p n -\(t,dq), 


( 2 . 21 ) 




{H\p*,Pi} - H\p i _ ll p*])f(t, pj-u dp*)f(t,p*,dpi) 
f(t, p i - 1 ,dp i ) 


Pn-i (t,dq) 


( 2 . 22 ) 


for i = 1,..., n — 1, and 




>n-i, P+]f{t, pn-i,dp+)p n -i(t, dq). 


(2.23) 


Remark. It is essential that the integral over B can be approximated by an integral 
over [x \,..., x n -i, L;po,..., p n - i, p+). The measure C*p below does not have any 
singular factors like 5(x n = L ). In particular, the result of integrating over B will 
partially cancel with the term arising from random particle entries. 


X 2 

▲ 



► X'l 


Figure 2. The deterministic flow <f> on the x-simplex is translation 
at constant velocity unless this translation crosses a boundary of 
A^. In the case n = 2 pictured above, points in Aq and Ai hit the 
boundary faces X\ = 0 and X\ = x’ 2 , respectively. The remaining 
portion of the simplex is mapped to the simplex minus the set B. 

The final ingredient for the proof of Theorem l2.1l is the time derivative of p,(t, dq ), 
which we now record. 
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Lemma 2.8. For any n > 0 and any 0 < s < t we have 

|| Hn(t, •) - fin(s, s)(C* Hn)(t , -)||tv = o(t - S ) (2.24) 

where the norm is total variation and (C* p n )(t, dq) is defined to be the signed kernel 


(C°£)(t,dp 0 ) | y (r/)(t,p,;-i,rfp.) 

£{t,dp 0 ) " f(t, Pi-1, dpi) 

n 

x e~ XL l A L{x!, ...,x n ) dx i • • • dx n £(t, dp 0 ) f(t, pj-i,dpj). (2.25) 

j =i 


The expression for the measure above is to be understood formally; the correct 
interpretation involves replacement, not division. All of the “divisors” above are 
present as factors of £(t, dpo) 11^=1 /(L Pj-i,dpj), and the fractions indicate that the 
appearance of the denominator in this portion is to be replaced with the indicated 
numerator. 

So that we know what to expect, before proceeding we note that when we sum 
over i in (12.251) . some of the terms arising from £° and £ K cancel. Namely, the 
bracketed portion of (12.251) expands as 


f H[p*,p 0 ]£(t,dp*)f(t,p*,dpo) 


-J H[po,P*]f{t,Po,dpi) 


£ 


£{t, dp 0 ) 

I{H[p*, Pi\ - g[pi-i,p*])/(<,/Ji-i,d/3*)/ft,p*,dpi) 

f(t, Pi-i, dpi) 

- J(H[pi,P*]f(t,Pi,dp*) - H[pi_i,p*]f(t, pi-i, dp*)) 


(2.26) 


The gain terms associated with the kinetic equations we leave as they are, but note 
that the “loss” terms telescope, and the above may be shortened to 


/ H[p*,po\£(t, dpj)f{t, p *, dpo ) 
i(t, dpi) 


- J H[p n ,p*]f(t,p n ,dp*) 


J{H[p*, P i] H\pi-i,p.])f(t, Pi-1, dp*)f(t,p*,d Pi ) 

(2.27) 


i— 1 


f(t, Pi-1, dpi) 


We are ready to prove our statistical characterization of the bounded system. 


Proof of Theorem \2.1l For times s and t with 0 < s < t < T, consider the difference 
f F{t, q) p(t, dq) - f F(s , q) p(s, dq): 


J l F (t, q) ~ F(s, g)] p{t, dq) + J F(t , q)[p{t, dq) - p(s, dq)] 

- J [F^, q) - F(s, g)] [p{t, dq) - p(s, dq)] = (I) + (II) + (III) (2.28) 

Since F{t,q) is uniformly continuous in t uniformly in q by Lemma 12.51 and p is 
a probability measure, (I) —> 0 as t — s —> 0. Using Lemma 12.81 and the fact that 
|F| < 1, (II) -A- 0 as t — s —> 0, and in fact 

/ Zf’ d7} / F(t, ,)(£•,)«.« (2.29) 
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as s —>• t— with t fixed. Using both Lemmas l2.5l and r2.81 we see (III) is o(t—s). Thus 
f F(t , q) pft, dq) is continuous in t; we will show additionally that it is differentiable 
from below in t with one-sided derivative equal to 0 for all 0 < t < T. In light of 
(12.291) . our task is to show that — f F(t,q) (£* p)it,dq) approximates (I) up to an 
oft — s) error. 

Using Lemma 12.61 we have the following approximation of the portion of (I) 
involving p n , with error bounded by C\ ,H'(P)[(t - s) 2 + {t- s)w{t - s)]: 


n— 1 


F(t, q)(l B (q) + 1 v(q)) - F(t, <j>\q) ^lcr(g) + Y Wg)J 
+ (t~ s)F(t, q) j H[p n , p+]f(t, p n , dp + ) 

-ft - s) J EF(t,e p+ <l>lq)H[p n , p + ]f(t, p n: dp+) 


(2.30) 


p n (t,dq). 


We have J[F(t, q)ly(q) — F(t, (j) t s q)lu(q)] p n ft, dq) = 0 because the deterministic 
How </>* maps U to V (modulo lower-dimensional sets) and preserve the Lebesgue 
measure in spatial coordinates. Making replacements using Lemma l2.7l and reorder¬ 
ing the terms, we find with an error bounded by CL, n ,\,H'(P)[ft — s) + wft — s)] 
that 


' F(t,q) - Fjs, q) 
t — s 


p„ft,dq) 


J F{t, e p+ q)H\p n -i, p+\fft, p n -\, dp + )p n -i{t, dq) 

- J E F(t, e p+ (j) T s q)H[p n ,p + }fit 1 p n , dp + )p n ft, dq) 

/77,-1 

r F(t q) ( ~ Pi-i, d P*) fit, p*,dpi 


(2.31) 


fit, Pi-i,dpi) 


H\p*,Po\£{t, dp*) fit, p*,dp 0 ) 


Pn-! (t,dq) 


lit, dpo) 

- J F(t,q) J H[p n ,p + ]f it, p n ,dp + )p n it,dq) 


for positive integers n. In the case n = 0, we have (p^q = q , and the approximation 
is 

J F it, g) _^i S,q ) p 0 it,dq) « j[Fft,e p „q) - Fft,q)]H[p 0 ,p*\poft,dq). (2.32) 

We observe that the bracketed portion of (12.311) nearly matches (12.271) , except that 
part involves p n -i and part involves p n . Furthermore, we have 


EF(t, t p+ <j> T s q)H[p n , p+\fit, p n , dp + )p n ft, dq) 

~ J Fit, e p+ q)H[p n , p+)fit, p n , dp + )p n it, dq) 


(2.33) 


< C H ’(P),n,Lit ~ S). 


















16 


KASPAR AND REZAKHANLOU 


This follows since r is conditionally independent of q given p n , and r E ( s,t ) a.s. 
so that all but 0(t — s ) volume of the x-simplex is simply translated by <j) T s to a 
region of identical volume. We make the replacement indicated by (12.3311 in (12.311) 
without changing the form of the error. 

For any positive integer N we define 

N 

r A r(t)= F(t,q) ^2 p n (t,dq). (2.34) 

^ n—0 


Summing (|2.31|) . (|2.32p , and our approximation for n n {t, •) — /i n (s,-) from (|2.27|) 
gives 


r N{t) ~ ^n{s) 


t — s 


N 


F(t,q) E 


i— 1 


{H[p*, pi] - H[p i ^ 1 , p*])fjt, p i - 1 ,dp*)fjt,p*,dp i ) 

f{t, Pi-i, dpi) 

H[p*,p 0 ]t(t, dp*) fit, /?*, dpo) 


(2.35) 


£(t, dp 0 ) J 

~ J F(t, £ P+ q)H[p N , p+\f(t, p N , dp+)p N {t, dq) 


p N (t,dq) 


with an o(t — s) error. Call the right side 7 jv(t), so that the left derivative of IV is 
7 n- Since T^r is a Lipschitz function, we deduce 

T n (T)~T n {0) = [ 7 N [t)dt. (2.36) 

J 0 


We have 7 ./v(t) bounded uniformly in t by 

, , .. 3H'(P)L N \ N+1 

IwWI S jjv _!), . 

and thus 7 jv —> 0 uniformly in t. It follows that 

r N (T)-r N (o)^o 

as N —> 00 . Recognizing the limits oiT N (T) and Tjv(O) as 

J G{q)p{T,dq) and J EG($%q)p(0, dq), 
respectively, we have verified (12.1411 and completed the proof. 


(2.37) 


(2.38) 

(2.39) 
□ 


Having shown the candidate measure p{t, dq) constructed using the solutions 
£{t,dpo) and f(t,p-,dp+) given by Theorem 11.61 our next task is to verify the 
latter, which we undertake in the next section. 


3. The kinetic and marginal equation 

The primary goal of the present section is to prove Theorem 11.61 concerning the 
kernel fit, p_, dp+) and marginal lit, dpo), so that the candidate measure pit, dq) in 
the previous section is well-defined and has the properties required for the argument 
there. 
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While we introduce our notation, we also discuss the intuitive meaning of the 
terms of our kinetic equation, comparing with that of Menon and Srinivasan [25] . 
Let us write C K for the operator on K. given by the right-hand side of (11.271) . 


(C K k)(p_,d P+ ) 


J{H[p*, p+] - H[p_, p*])k(p_,dp*)k(p*,dp+) 

J H[p+,p*]k(p+, dp*) ^ J H [P-,P*]k(p~, dp*) 


k(p-,dp+), 

(3.1) 


so that the kinetic equation is ft = C K f. 

The first term, which we call the gain term, corresponds to the production of a 
shock connecting states p _ and p+ by means of collision of shocks connecting states 
P- , p * and p*, p+. Such shocks have relative velocity given by 


H[p*,P+\ ~ H[P-,P*\ = H[P-,P*,P+](P+ ~ P-)- (3-2) 

In the Burgers case, the second divided difference is constant, and the 

above is proportional to the sum of the increment p+ — p_, analogous to mass in 
the Smoluchowski equation. 

The second line of HO we call the “loss” term, though this need not be of definite 
sign. To better understand this, note that when we have proven Theorem 11.61 we 
will know that /(£, p_,dp+) has total mass which is constant in (f, p_). In this case 
the loss term may be rewritten equivalently as 


J(H[P+,P*\~ H[p-,p+])f(t, p+, dp*) 

- J(H[P-,P*\ - H[p-,p+))f(t, p -, dp*) 


f(t,P-,dp+), (3.3) 


which corresponds precisely with the kinetic equation of [25] . The meaning of 
the first line of (G3D is clear: we lose a shock connecting states p-,p+ when a 
shock connecting p+,p* collides with this, and the relative velocity is precisely 
H[P+,P*\ - H[p-,p+]. 

The second line of (13.31) is less easily understood. One would expect to find 
here a loss related to a shock connecting p_, p+ colliding with p* , p_ for p* < p_, 
particularly if we were viewing / as a jump density as [25] views its corresponding 
n(t, y, dz). Viewed as a rate kernel , it is not clear that / should suffer such a loss: 
in some sense we must condition on being in state p_ for /(t,p_, •) to be relevant 
at all. At this point we cannot offer an intuitive kinetic reason for the second line 
of (031), but in light of the rigorous results of this article we can be assured that 
this is correct for our model. 


Remark. The form (EU) seems preferable in the more generic setting, since—as 
pointed out by the referee—the resulting equation f t = C K f has the property that 

Ht,P-)= J f(t,P-,dp+) (3.4) 

is (formally) conserved in time for each p_, without assuming that A(0 ,p~) is con¬ 
stant. This observation will be important in attempts to generalize the results of 
the present paper. 

We return to the task at hand, showing existence and uniqueness of /(£, p_, dp+). 
The argument proceeds through an approximation scheme for exp(A H'(P)t)f(t, •, •), 
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where we can easily maintain positivity. To avoid endowing /C with a weak topology, 
we show directly the approximations are Cauchy, rather than appealing to Arzela- 
Ascoli. 


Proof of Theorem \1. (A part I. Write c = XH'(P) and consider for each positive 
integer n the continuous paths h n : [0, oo) —> K, defined by h n ( 0) = g for all n and 

h n (t) = e~ cj/n (C K h n ){j/n) + ch n (j/n ), t e (j/n, {j + 1 )/n). (3.5) 

We claim the following properties of h n : 

(a) h n (t) > 0 for all t > 0, 

(b) for each t > 0 the total integral f h n (t, p-,dp + ) is constant in e [0, P), 

(c) for all t > 0 

\\h n (t)\\ < A(1 + c/n) rntl < Ae c(t+1) =: M(t), (3.6) 


and 

(d) h n (t,p-, [0,p_] U { P }) = 0 for all t and all p- € [0,P). 

Since h n is piecewise-linear in t and the properties above are preserved by convex 
combinations, it suffices to verify this at t = j/n for all integers j > 0. 

We proceed by induction. Abbreviate hf = h n (j/n). The j = 0 case holds by 
the hypotheses on g. Assume now that the claim holds for j, and consider j + 1. 
We have 


n(h ™ + 1 - h]) = e~ cj/n / {H[p*,p + ] - H[p_, p*])ft.™(p_, dp*)h]{p*, dp+) 

J H ^ p -' p ^ h ^ p -' dp *^) h ?( p ~’ dp +) 

J tf[p+,P.]fr"(p+,dp.)^ h](p-,dp+) 

+ chj(p-,dp+). 

In particular, for any p + we have 

c-e- c V n j H[p + ,p*]h™(p + ,dp*)>c-e- c i/ n H'(P)\\h]\\ 

> c - e~ cj / n H'{P) A(1 + c/n) j 


(3.7) 


(3.8) 


= c 


l-e~ cj/n {l + c / n y 


> 0 , 


using property (c) for case j; thus n(/i" +1 — hj) > 0. This and (a) h™ > 0 give 
Pyfi > 0, property (a) for case j + 1. Furthermore, the total integral of the right- 
hand side of m is exactly 


: J h](p_,dp + ) =c\\h r - 


(3.9) 


since (b) implies the bracketed portion integrates to 0. Also using (b) for case j, 
the change in the total integral from K'f to /t" +1 is independent of p_, verifying (b) 
for case j + 1. Using (c) for case j we find 

K+rII < \\h?\\+ cn- 1 ||/i”|| = (l+c/n)\\h]\\ < A(1 +c/ n y, (3.10) 

so that (c) holds for case j + 1. Finally, we observe that property (d) for hf implies 

C K h](p-,[0,p.\U{P}) = 0, (3.11) 
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and thus (d) holds for case j + 1. 

We now consider the matter of convergence. Pairing C K with any k\,k 2 € K. and 
measurable J(p+ ) with | J| < 1, we find that 

II C K h - 1| < 3H'(P)(\\ki\\ + ||fe||)||fe - fc 2 ||. (3.12) 

In particular C K is Lipschitz on bounded sets in /C. For brevity write 

C h k(s, p-,dp + ) = e~ cs C K, k(p- 1 dp + ) + ck(p-,dp + ). (3.13) 

We compute for any m < n and t > 0 

|| h m (t) - h n (t )|| < [ || h m (s) - C h h m (s)\\ + \\C h h m (s) - C h h n (s)\\ 

Jo 

+ \\C h h n (s) — /i"(s)|| ds. (3.14) 

Observe that 

||/i n (s) - C h h n (s)\\ = \\C h h n (n~ 1 |yisj) - C h h n (s )|| 

< Ue-cn-'L^J^^^-i^sj) -e~ cs C K h n {s)\\ 

+ c||/i n (n _1 |_nsj) - h n (s) || 

where 

||e _cn_1 L" s J£ K /i"(n _:L LnsJ) — e~ cs C K h n (s)\\ 

< e~ cs ||£ K /i"(n _1 |_nsJ) — C K h n (s)\\ 

+ (e- cn_1 L ns J - e~ cs )\\C K h n (n- 1 |nsj)|| 

< e- cs 6H\P)M{s)\\h n {n- 1 [ns\) - /i"(s)|| 

+ e~ cs e cn ~ 1 n~ 1 3H' (P)M(s) 2 

< 6ce c ||/i n (n _1 |_ns_|) — /i"(s)|| + 3cAe C5+3c n _1 , 

We have also 

||/i n (n^ 1 (nsj) — /i n (s)|| < n _1 ||£ ,l /i r! '(n _1 LnsJ)|| 

< n~ 1 (e~ cs 3H , (P)M (s) 2 + M(s)) 

= n~ 1 (3Xce cs+2c + Ae c(s+1) ). 

Combining these, there is a constant C\ so that 

\\h n (s) ^ C h h n (s)\\ < C\n~ 1 e cs 

for all positive integers n and times s > 0. 

The remaining term in the integrand of (13.141) is 

|| C h h m (s) - C h h n (s) || < e -C5 ||£ K /i m (s) - C K h n (s)\\ + c||/i m (s) - h n (s)\\ 

< (e~ cs 6H'(P)M(s) + c)\\h m (s) - h n (s)\\ (3.19) 

= (6e c + c)\\h m (s) - /i"(s)|| = C 2 \\h m (s) - h n (s)\\. 

All together, 

\\h m (t)-h n (t)\\< [ C 2 \\h m (s)-h n (s)\\+C 1 (n- 1 + m~ 1 )e cs ds, (3.20) 

Jo 

and by Gronwall 

|| h m (t) - h n (t) || < CAn~ l + m~ 1 )e C2t ^——^. (3.21) 

c 


(3.15) 

(3.16) 

(3.17) 

(3.18) 
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From this we see that h n is Cauchy, hence convergent to a continuous h : [0, oo) —» K, 
satisfying 

h(t) = g + f [e~ cs (C K h)(s) + ch(s)]ds (3.22) 

Jo 

for all t > 0, having properties (a,b,d) above and ||/i(i)|| < M(t). 

We define f(t) = e~ ct h(t), finding that 

• t i-A f(t) £ JC + is continuous; 

• f H C K f(t) £ K. is continuous (using the local Lipschitz property of £ K ), 
and hence Bochner-integrable; and 

• / = C K f with /(0) = g by Leibniz. 

Using the local Lipschitz property of C K , the solution is unique. Properties (a,b,d) 
above hold for /, and (b) in particular gives for each fixed p_ 

j t J f{t, dp+) = J (C K f)(t, p_, dp + ) = 0. (3.23) 

Thus f f(t, p-,dp + ) = A for all f > 0, 0 < p_ < P. □ 

We now turn to £(t, dpo), which is intended to serve as the marginal at x = 0 for 
our solution p(x , t) for fixed t > 0. We define a time-dependent family of operators 
C° acting on measures v{dpo) in M. by 


(C°v)(t,dp 0 ) = / H[p t ,pQ]v(dp*)f(t,p t ,dp 0 ) 


H[p 0 ,P*\f(t, Po,dp*) 


v(dp 0 ). (3.24) 


Again the integration is over p* only, and for each t we have •) £ A4. Our 

evolution equation for £(t, dp o) is the linear, time-inhomogeneous £ t = C°l. 


Proof of Theorem 1 1. 61 part II. Note that (13.2411 is exactly the forward equation for 
a pure-jump Markov process evolving according to a time-varying rate kernel 

H[Po,P+]f(t,po,dp+), (3.25) 

so we could obtain existence and uniqueness of solutions along these lines. We 
outline an argument similar to that employed for /, for the sake of completeness. 
Define for each n the continuous path h n : [0, oo) —» A4 with h n ( 0) = <5o and 

h n (t) = e~ cj/n C°h n {j/n) + ch n (j/n), t £ {j/n, (j + l)/n). (3.26) 

Proceeding as in the proof for / we verify that for each j the hf = h n (j/n) are 
nonnegative and, since C°h n (j/n ) has total integral zero, ||/i"|| < (1 + c/n)^ rat l. 
Observing that C {) is linear and bounded uniformly over t, 

sup \\{C°v)(t, - JUtv < 2XH'(P)\\u\\ tv , (3.27) 

t 

we easily show h n is Cauchy on bounded time intervals, with limit h : [0, oo) —>• ]£[+ 
satisfying 

h(t) = Sq + [ [e~ cs (C°h)(s) + ch(s)]ds. (3.28) 

Jo 

Take £{t) = e~ ct h(t) to find that £(t) > 0 solves 

£{t)=5 0 + [ ( C°£)(s)ds . (3.29) 

Jo 
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Uniqueness follows using (13.271) . Finally, since C Q v as zero total integral for any v , 
we find that f £(t, dpo ) = 1 for all t. □ 

We close this section with a remark concerning the kinetic equation f t = C K f 
without the assumption that the initial kernel g has bounded support [0, P], When 
the initial condition £(x) is unbounded, growing for example linearly in the case 
of quadratic iJ, the solution to the scalar conservation law p t = H(p x ) on the 
semi-infinite domain should blow up in finite time. We likewise expect that the 
solution / to the kinetic equation will blow up, but control of certain moments 
prior to this blow up may allow us to run the argument of Theorem 12.11 with only 
superficial changes. At the moment we lack the sort of estimates one obtains in the 
Smoluchowski case (e.g. coming from a closed equation for the second moment), 
but we hope to revisit this in future work. 

4. Conclusion 

We review what has been accomplished: using an exact propagation of chaos 
calculation for a bounded system with a suitably selected random boundary condi¬ 
tion, we have derived a complete description of the law of the solution p{x , t) to the 
scalar conservation law p t = H(p) x for x,t > 0 when p(x, 0) = £(a;) is a bounded, 
monotone, pure-jump Markov process with constant jump rate. Notably, we have 
recovered the Markov property of the solution and a statistical description of the 
shocks simultaneously. This may be regarded as both a strength of our present 
analysis and a shortcoming of our present understanding (a soft argument for the 
preservation of the Markov property in the particle system would be illuminating). 

We emphasize how our approach using random dynamics on a bounded interval 
has made things considerably easier: by constructing a bounded system which has 
exactly the right law, we are relieved of the burden of determining precisely how 
wrong the law would be with deterministic boundary. 

Our result lends additional support to the conjecture of Menon and Srinivasan 
m and we hope that the sticky particle methods described herein will provide 
another approach, quite different from that of Bertoin [4], which might be adapted 
to resolve the full conjecture. One of the authors is attempting a similar particle 
approach in the non-monotone setting, and hopes to report on this in future work. 

Acknowledgments. The authors thank the anonymous referee for suggestions 
improving our presentation and for pointing out a recent reference. 

Appendix A. Proofs of lemmas 

Proof of Lemma I A 51 Boundedness of G is obvious from its definition as a Laplace 
functional. For continuity, we begin by considering G(cf>Qq) with q fixed. As t varies, 

q varies continuously, with p -values unchanging and x-values changing with rate 
bounded by H'(P ), except for finitely many times when collisions occur. Note 
G{4> t 0 q) depends on <f>Qq via 7r(</>gg, •), and the former varies continuously across these 
collision times. More concretely, because J G C([0, L]), J is uniformly continuous, 
and so 


wj(S) = sup{| J(a;) - J(y)\ : x, y G [0, L\, \x - y\ < <5} 


(A.l) 
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has the property that wj(S) —> 0 as S —>• 0+. Using the above, 

n 

llogG&lq) -logG(q)\ < ~ Pi-i)wj(tH'(P)) < Pw(tH'(P)). (A.2) 

t=1 

Since the exponential is uniformly continuous for nonpositive arguments, we find 
that G{<f>\q) is continuous at t = 0 uniformly in q. Write 

wg(S) = sup{|G(0o<z) - G{q)\ : q £ Q, 0 < t < <5}, (A.3) 

and observe that wg(S) —> 0 and S 0+. 

Suppose now that 0<s<f<T;we compare EG(4>;fg) and E G(<j>fq). Write 
9 = t — s. We have G($Jg) = G(d>^_ e $J’ _6 'q) a.s., coupling the random entry 
processes on matching intervals. Then 

E| G(^_ e ^- e q) - G((fc_ g *T- 9 q) | < C XiH ' {P) 0, (A.4) 

because the random entries occur at a rate bounded by XH'(P) and |G| < 1. Using 
the time homogeneity and continuity for the deterministic flow, 

E| G{<fe_ e $ T s - e q) - G(^~ 0 q)\ < w G {0). (A.5) 

So, at the cost of an error bounded by C\^h’{p)9 + wg(9 ), we compare E G{<&Jq) 
with EG($J~ 6l g) instead. Now the configuration q is flowed on time intervals of 
equal length, but with random entry rates that have been shifted in time. 

The random entry rate is given by H[p n , p+\f(r, p n , dp+) and / is TV-continuous 
in r; we have 


II f(r,p n ,dp+) - f(r - 0,p n ,dp+)\\ 


t r—0 


(G K /)(r, p n ,dp + )dr 




e. 


Let us define three kernels for r £ [0, T — t] according to 

f(r, P-, dp+) = f(s + r, p-, dp+) A f(t + r, p -, dp+) 
f s ( r i P-i dp+) = f(s + r,p-,dp + ) - f(r,p-,dp+) 
f t { r i P—i dp+) = f{t + r,p-,dp+) - f(r,p-,dp+), 


(A.6) 


(A.7) 


where the minimum (A) of two measures is defined as usual by choosing a third mea¬ 
sure with which the former two are absolutely continuous, and taking the pointwise 
minimum of their Radon-Nikodym derivatives. That this extends measurably to 
the parametric case (i.e. involving kernels) is immediate from a parametric version 
of the Radon-Nikodym theorem [2SJ Theorem 2.3]. 

Using the above we construct a coupled random entry process for times r £ 
[0, T—t], namely let {( s ,( t )( r ) be the pure-jump Markov process started at (p n , p n ) 
from q and evolving according to the generator C z acting on bounded measurable 
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functions J(y,z ) defined on [0, P] 2 according to 

{C z J)(r,y,z)= 1 (y = z) J{J{p + ,p + ) - J(y,z)}H[y,p + ]f(r,y,dp + ) 

+1 (y = z) j{J(p+,z) - J{y,z)}H[y,p + ]f s (r,y,dp + ) 

+1 (y = z ) j{ J{V,P +) - J{v,z)}H[z,p + ]f t (r,z,dp + ) 

+1 {y ± z ) JU(p+,z) - J(y, z)}H[y, p+\(f + f s )(r,y, dp+) 

+ Mv ¥= z) f {J(y,P+) - J(y,z)}H[z,p+\(f + f t ){r,z,dp+). 

(A.8) 

Taking J(y , z) which does not depend on z we find that ( s (r ) for re [0, T — t] has 
the same law as the random boundary for and likewise taking J which does 

not depend on y we find the £ t (r) has the same law as the random boundary for 
<1>J g, verifying the coupling. 

On the diagonal y = z, the rate at which C z causes jumps that leave the diagonal 
(the second and third lines of (|A.8|I ) is bounded by 


Ca ,h'{P)0i 


(A.9) 


and the probability that such a transition never occurs in a time interval of length 
T — t is bounded below by 


exp [—C\ tH '(p){T — t)Q] > exp[—Cx, H '( P )Td\. 


(A.10) 


So, coupling the random entry dynamics, we find that 9 q = <f 'Jq with proba¬ 
bility at least exp[— Cx,h'(p)T0\. Putting the above pieces together, we find 

\EG(^q)~EG(^q)\ 

< C x ,H>(p){t - s) + w G (t - s) + (1 - exp(Cx,H'(p)T{t - s))), (A.ll) 
and the proof is complete. □ 


Proof of Lemma \2.6\ Using the Markov property of the random flow $, we have a 
functional identity: 

F(s,q) =EG($fa>‘g) = E F(t,& s q). (A.12) 

Let q be fixed, and consider the following events, whose union is of full measure for 
computing the final expectation above: 

£o = {no entry at x = L in (s, t)} 

£\ = {at least one entry at x = L in (s, t)} 

Observe that on £q we see only the deterministic flow <j) over the time interval (s, t ): 

F(t,& a q)leo=F(t,<l>*q)leo, (A.14) 

and this occurs with probability 


(A.13) 
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We prefer an expression evaluating / at only a single time, and Taylor expand the 
exponential around zero. 


J H[Pn,P+]f(r, p n ,dp+) dr- (t-s) j H[p n , p+\f(t, p n , dp+) 


III H [Pn,P+](F- K f)(r',p n ,dp + )dr' dr 
J S Jt J 


< C\,H'(P){t — s) 2 , (A.16) 


so that 


P(£o) = 1 -(t-s) j H[p n , p+\f(t, Pn , dp+) + 0((t - s) 2 ) (A.17) 


and by exhaustion 

P(£i) = (t-s) J H[p n , p+\f(t, p n , dp+) + 0((t - s ) 2 ), (A.18) 

with both errors bounded uniformly over Q. On £\ write r for the first time a 
random entry occurs for $*, and p + for the new boundary value, noting that the 
distribution of r depends only on q through p n and that the law of is determined 
by /(t, p n , ■)■ We have r £ ( s,t ) so that 


F(t,& a q)l £l = F(t,& T e p+ 4> T s q)le x . (A.19) 


Using the strong Markov property for the random boundary at the stopping time 

t, 

EF(t,^ t T e p+ (l) T s q)l £l =EF(T,e p+ (j) T s q)l£ 1 . (A.20) 

Since P(£i) = 0(t — s ), we can afford to make o(l) modifications to this. Write 
w(S) for the modulus of continuity of F(t, q) in time, according to Lemma l2.5l Now 
t G (s, t) a.s. on £i, so 


I ^F(T,e p+ (/) T s q)l £l -EF(t,e p+ 4>lq)l £l \ < P(£i )w(t-s). (A.21) 


Next we modify the distribution from which p + is selected; at present, p+ is 


selected according to the random measure 

H\p n ,P+]f(r, p n ,dp+) 


f H[p n ,P+\f(T, Pn,dp+)' 

Since r G (s,t) a.s. on the total variation difference 

\\f(T,p n ,dp+) - f(t,p n ,dp+)\\ < Cx,H'(P)(t - s) a.S. 
Let us write p + for an independent random variable distributed as 


(A.22) 
(A.23) 


H{pn,p+)f(t,p n ,dp + ), (A.24) 

(note t instead of r), normalized to have unit mass. Then 

\EF(t,e p+ <plq)l £l — EF(t, e p+ (j> T s q)l £l \ <C x , H >( P )(t-s) P(fi). (A.25) 

Note that P(£i) cancels with the normalization in (IA.22D . up to an 0(t — s) error, 
and that is independent of £\. With error bounded uniformly over q G Q by 
C\,H’(p)[(t — s) 2 + (t — s)w(t — s)], we have 


EF(t, $‘g) « F(t, <j>\q) ^1 -(t-s) J H[p n , p+]f(t, p„, dp+) 


+ (t — s) 


K[F(t,e p+ cl) T s q) | £i]H[p n , p + \f(t, p n ,dp + ), (A.26) 
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where the remaining expectation is now taken only over r conditioned on £\, which 
is therefore distributed over [s,f], □ 

Proof of Lemma UfTA On any of i = 0,..., n — 1, we have a = a(q) < t. For 
such q we have 


I F(t,4>U) ~ F (t, <PsQ)\ 

< I ~ - F (t,<l>s<l)\- ( A - 27 ) 

The first absolute difference is the error induced by conditioning on zero particle 
entries on the interval [s, cr], which costs less than A H'(P)(t — s). Lemma l!T5l bounds 
the second of these by w(t — s). Using (IA.271 1 we can replace the integrand F(t, 4>\q) 
over Ai with F(t, 4>fq), i.e. the projection of q onto (9jA£) x A^ +1 in the direction 
of the velocity. 

Write Vi = —H[pi-i, pf\ for the velocities. We can integrate over Aq using as a 
parametrization for the x-variables 


(0,x 2 ,..., x n ) (~0V1,X 2 -0v 2 ,...,x n - 0v n ), (A.28) 


which has absolute Jacobian determinant —vi = H[po,pi\. Namely, 


F {t,<jf s q)p, n {t,dq) 

1 

n f'(t — s)A\(L — Xn)/Vn\ 


I A^_ 1 xA ^ +1 Jo 


F(t,(0,x 2 , . .. ,Xn, po, ■ ■ ■ , Pn)) 


H[p 0 , pi] dd dx 2 • • • dx n £(t, dp 0 ) f(t, pi-i,dpi). 

(A.29) 

On all but a portion of A 0 with x-volume bounded by — s) 2 we have 

t — s < |(L — x n )/v n \, and we can extend the upper limit of integration to 9 = t — s 
with and error bounded by this. Noting that the configurations 


(0, x 2 ,...,x n ;po,...,p n ) and (x 2 ,..., x n \ pi,..., p„) (A.30) 


are equivalent under our sticky particle dynamics and relabeling, we obtain our 
approximation for the integral over Aq. 

The argument for Ai, i = 1,..., n — 1, is similar. We use the mapping 


(xi ,...,xi,x i+2 ,...,x n ,6) i y (xi -Qv i,...,Xi - 6v it Xi - 0v i+1 ,... ,x n ~0v n ), 

(A.31) 

with absolute Jacobian determinant Xj—Wj+i = H[pi , pi+i]— H\pi_i, pf\, and replace 
the upper limit of integration for 6 with t — s at the cost of an 0((t — s) 2 ) error. 
Relabeling gives the indicated approximation. 

The range of 

<f> s s +6 (x l, • • •, x n -!,L- po,..., p n ) (A.32) 

where (xi,... ,ar„_i) G A£_ 1; (p 0 ,... ,p n ) G A^ +1 , and 


0 <0 <{t 


x x i 

s A —A 

\vi\ 


x 2 - X\ 
\v 2 - Vi| 


A • • ■ A 


\Vn 


%n —1 
Vn—l | 


(A.33) 


is B, modulo sets of lower x-dimension. Using Lemma f2.5l we can replace F{t,q) 

by 

(A.34) 


F{t, (Xl, . . . , Xji — 1 , L’, po, , p n )) 
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with pointwise error over B bounded by C\ t H'(p) \w{t — s) + (t — s)]. For all but an 
0{(t — s) 2 ) portion of the x-volume, the value t — s is the minimum above, and as 
before we obtain the indicated approximation of the integral over B. □ 


Proof of Lemma 1 2. #1 Essentially we are verifying the Leibniz rule, but we are un¬ 
able to find a version of this to cite for kernels. We first obtain quantitative control 
over our linear approximations of / and i. Namely, fix any measurable | J| < 1. We 
have 

J J(p+)[f(t,p~,dp+) - f(s,p-,dp+) - (f- s)(C K f){t,p-,dp+)} 

= [ J(P+)[(£ K f)(T,P-,dp+) - (C K f)(t,p-,dp+)]dT. (A.35) 

J S 

The difference of £ K / at different times can be expressed in terms of / and H again, 


(£ K f)(r, p-,dp+) - (C K f)(t, p-,dp+) 

= j(H\p*,P+\~ H\P-,P *]) Q" {£ K f)(T l ,P-,dp*)dT , ' S j f(r,p*,dp+) 
+ J(H[P*,P+] ~ H[P-,P*])f(t,P-,dp*) {C K f)(r', p*, dp + ) dr' 

- J J (H[P+, P*](£ k )(t', p+, dp*) 


(A.36) 


- H[P~, P*\{£ K ){t' , P-,dp*) ) dr 1 

J H[P+,P*]f(P P+, dp*) 


/( t , p-,dp*) 


- H[P-,P*\f( T ,P-,dp *) 


j\c K f)(T', P -,dp + ) 


Noting that ||£ K /|| < 3A 2 H'(P), we integrate over p +1 then p*, and find that (IA.35D 
is bounded by 


9H'(P) 2 X 3 (t — s) 2 . 


(A.37) 


Next, for any \ J\ < 1, 


J{Po)[t(t,dp 0 ) -£(s,dp 0 ) - (t - s){C°e)(t,dp 0 )} 

= [ J (Po)[(£°^)(T,rfpo) - (C°t)(t,dpo)]dT, (A.38) 

J S 
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and 


(£°f)(T, dp 0 ) - (. C°£)(t, dp 0 ) 

= j H[p*,po\ (C°£)(t' ,dp*)dr^J f(T,p*,dpo) 

+ J H[p*,p 0 ]£(t,dp„) (C K f)(T',p*,dp 0 ) 

H[p 0 , p*\{C K f){r', p 0 ,dp*)\ dr' 


(A.39) 


LJt 


£(r,dp 0 ) 


H[po,P*]f(t,Po,dp*) 


(£°£){t', dp 0 ) dr'. 


J Jt 


We have ||£°£|| < 2A H'(P), so by integrating over p 0 and then p* we find (IA.38I) is 
bounded by 

5H'(P) 2 X 2 (t — s) 2 . (A.40) 

Returning to the problem of establishing our Leibniz rule, note that e~ AL l ^(dx) 
factors from both p n and C* p n . It will therefore suffice to obtain a bound on the 
(po, ■ ■ ■, Pn) portion. 

We argue by induction. In the case n = 0, we have only the difference in (1A.38|) . 
and the result holds. Now suppose that the result holds for case n , and consider 
n + 1. Choose as a test function J(po, ■ ■ ■, Pn+i), which is measurable and has 
| J| < 1, and integrate it against 


n+1 


n+1 


£{t, dpo) / (t, Pj—i, dpj) - £(s,dp 0 ) f(s, pj-!, dpj) 


7 = 1 


P n 

= / £{t,dp 0 ) Y[f(t,pj-i,dp j )[f(t,p n ,dp n+1 )-f(s,p n ,dp n+1 )\ 
3=1 


+ 


t(t,dp 0 )Y[f{t,Pj-i,dp j ) — £(s, dpo) f(s, pj-i,dpj) 
7=1 7 = 1 


/(fi) Pm dp n + 1). 

(A.41) 

For each po,...,p ra , the function J(po, ■ ■ ■, Pn+i) is bounded and measurable in 
Pn+ 1 , so f(t , p n , dp n+ 1 ) - f(s, p n , dp n+ 1 ) can be replaced with 

(t~s) J J(po,...,p n+1 )(C K f){t,p n ,dp n+1 ) (A.42) 


plus an error no larger than 9H'(P) 2 X 3 (t — s) 2 , which when integrated over the 
remaining variables grows by a factor A™. 

In the third line of (IA.41I) . noting f J(po,..., p n +i)f(s, p n , dp n +\) is bounded 
and measurable in (po, ..., p n ), we apply the inductive hypothesis, replacing the 
bracketed difference by 


(f-s) 


(■ C°£)(t,dpo ) (£+)(£ p^+p,) 

£{t,dp 0 ) “ f(t, pi-i, dpi) 


£(t, dpo) n /(£ Pj-i)dpj) 
7=1 


(A.43) 


plus an o{t — s) error. After doing this, again noting J(po,..., p n + i) is measurable 
in p n+ 1 for each fixed po,---, Pm we replace f(s, p n , dp n+1 ) with f(t, p n , dp n+ 1 ) at 
a cost of 3A 2 H’(P)(t — s), which gets multiplied by the other factor (t — s). 
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Adding the modified versions of these two lines of (IA.41|) . we find exactly the 

po ,..., p n +1 portion of C*p n +1 plus an o(t — s) error. □ 
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